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Concept of transverse deflection probability of a parton that travels through strongly interacting 
medium, recently introduced by D’Eramo, Liu and Rajagopal, have been used to derive high energy 
evolution equation for the jet quenching parameter in stochastic multiple scatterings regime. Jet 
quenching parameter, q{x), appears to evolve with x, with an exponent 0.9tis, which is slightly less 
than that of xQ{x) where G{x) is the gluon distribution function. 
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I. INTRODUCTION 

The heavy-ion collision programs at CERN’s Large 
Hadron Collider have opened new access for exploration 
of extreme hot and dense nuclear matter. Precision to¬ 
mography of the nuclear matter are now become feasible 
as accurate test of the underlying quantum chromody¬ 
namic (QCD) theory. This is instrumental in discov¬ 
ering yet unexplored characteristics of various nuclear 
effects and collective phenomena that the nuclear mat¬ 
ter may possess. One possibility is to explore the QCD 
scale/energy evolution of various observables in this ex¬ 
treme ambiance. Advancement of study for hard sector 
observables at the LHC elevated the medium modifica¬ 
tion of high energy jets as prevailing topic of investiga¬ 
tion. In this context study on scale/energy evolution of 
the jet quenching parameter, attributed as the stopping 
power of the medium for certain probe of the medium, 
is now viable. High energy quarks and gluons passing 
through the interacting nuclear matter have their trans¬ 
verse momentum distribution broadened due to multiple 
scatterings with the constituents of the medium. While 
travelling through the strongly coupled medium the hard 
parton loses energy as well as its direction of momen¬ 
tum changes. Change in the direction of momentum is 
referred as ‘transverse momentum broadening’ for the 
travelling parton. In the context of jet-medium inter¬ 
action the broadening refers to the effect on the jet when 
the direction of the momenta of an ensemble of partons 
changes due to the random kicks. Even though there is 
no apparent change in mean momenta, the spread of the 
momentum distribution of individual parton within that 
ensemble broaden. 

Evolution of the momentum broadening was first stud¬ 
ied by Liou, Mueller and Wu Q by introducing radiative 
modification over the leading order momentum broaden¬ 
ing effect. The authors shown that average momentum 
broadening (p^) has both double and single logarithmic 
terms. Both the double logarithmic terms, ln^(L/^o), and 
single logarithmic terms, ln(L/(o) are coming from gluon 
radiation induced by the medium interactions. Here 
length of the nuclear matter is L and Iq being the size 
of constituents of the matter. Their estimation shown 
that the radiative contribution to be a sizeable correc¬ 


tion to the nonradiative leading value of {p^}. Later an 
evolution equation have been obtained for the inclusive 
one-gluon distribution, through the concept of classical 
branching process and cascade of partons This ex¬ 
plicitly takes into account the dependence of the observed 
gluon spectrum upon the energy and the transverse mo¬ 
mentum. The explicit transverse momentum dependence 
of the splitting kernel then enable one to identify large 
corrections to the jet quenching parameter. Subsequent 
studies on non-linear evolution leads to prescribe the 
renormalization of the jet-quenching parameter [^, Q. 

In this paper, in order to study the energy evolution 
of jet quenching parameter, we have adopted the idea of 
transverse deflection probability of a parton, that travels 
through the nuclear medium. Following a recent work 
by D’Eramo, Liu and Rajagopal [I| we then relate the 
momentum broadening to the S'-matrx of the nuclear in¬ 
teraction for a dipole. Balitsky-Kovchegov equation as 
the evolution equation of the S'-matrix is then used to 
derive high energy evolution equation for the jet quench¬ 
ing parameter in stochastic multiple scatterings regime. 
The known result of double log enhancement emerges 
as a special case in the limit when the single scatter¬ 
ing is only contributing. Power-counting techniques bor¬ 
rowed from Soft-Collinear-Effective-Theory (SCET) (Tfll - 
[ 2 ^ have been used to identify the leading contributions in 
the stochastic multiple scatterings region. For an almost 
constant q{uj) we recovered the double log result (in the 
limit t qL) first derived in Q and subsequent other 
studies [3-01 • We have also shown that the double log en¬ 
hancement diluted when we go beyond single scattering 
limit as previously argued in Q. Jet quenching param¬ 
eter q{x) is found to evolve with x, slightly weaker than 
xQ{x) where G(x) is the gluon distribution function. 

II. PROBABILITY DENSITY FOR 

TRANSVERSE DEFLECTION P(k±) 

The transverse momentum broadening of a parton can 
be studied by introducing concept of a probability den¬ 
sity, denoted in this article as P(k±). It signifies the 
probabilistic weight for the event where, after travelling 
a medium of length L, amount of transverse momentum 
that the parton acquired is k± [l|. The probability den- 
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sity P{k±) is chosen to be normalized in the following 
way, 

Using this probability density one can quickly estimate 
mean transverse momentum picked up by the hard par- 
ton per unit distance travelled, 


q = 


(kl) 


1 


47rL 


J dki ki p{k±). 


( 2 ) 


This defines the jet broadening (or quenching) parameter 
q that may have intrinsic dependence on x and Q of the 
probe through P{k±). 

Its important to make certain caveat when defining 
the jet quenching parameter, q, as done in Eq. ([2|). In 
absence of any finite upper bound to the integral, the 
right-hand side of Eq. ([2]) could be diverging, e.g, when 
P{k±) behaves as power law, P{k±) ^ 1/fc^ in the event 
of lowest order perturbative gluon production at large- 
k±. This gives a logarithmic UV divergence, resulting 
both average transverse momentum (fc^) as well as jet 
quenching parameter, q, infinite. Although, in case of 
multiple scattering with Sudakov like form factor or in 
the event of multiple stochastic scatterings where P{k±) 
takes the form of an exponentially damping factor, right 
hand side of Eq. m should be finite even without any 
apparent upper bound in the integral. 

In this paper we will discuss the event of multiple 
stochastic scatterings, where the probe receives random 
transverse kicks and P{k±) takes the form of a Gaussian 
with a variance of qL/2, 


47r f k'^\ 

P(k^) = ^ exp (-- 


( 3 ) 


In this study we are looking for transverse momen¬ 
tum broadening of the hard parton that have initial light 
cone momentum p{p'^,p~,p±) ^ (5(0,1,0) and enters 
in a brick of strongly interacting medium of length L. 
In order to do the power counting, we have introduced 
the dimensionless small parameter A. Power corrections 
in A to some hard process are generally suppressed in 
the presence of a hard scale, » ^qcd^ which act as 
base for the power corrections. This is a concept bor¬ 
rowed from soft collinear effective theory (SCET) studies 
In SCET studies the power counting protocol 
is as follows, terms that are sub-leading in two order 
of magnitude can be dropped but the terms that are 
suppressed by one order of magnitude should be kept. 
While travelling through nuclear medium the hard par- 
tons interacts repeatedly with the Glauber gluons which 
scale as (5(A^,A^,A). After first a few scatterings the 
hard parton’s momentum becomes of order (5(A^,1,A), 
though still quite collinear. In this scenario of high en¬ 
ergy regime where Glauber gluons are mostly effective, it 
has been shown in Ref. [I| that P{k±) can be expressed 
as the Fourier transform in r± of the expectation value 


of two light-like path-ordered Wilson lines transversely 
separated by rj_. 


where. 


P{kx) = J (frx S{rx) 


5(r^) = —(Tr[W(0,r^)W(0,0)] ) , 


( 4 ) 


( 5 ) 


with, 

W(y+,?/±) = P < exp 


ig / dy A+ {y+,y ,y±) 
Jo 


( 6 ) 


Within a dipole picture transverse separation yj_ — y'j_ 
can be taken as the transverse size of the dipole rj_. The 
length of the medium is L = L~ 1^/2. In this work we 
have assumed that P{k±) as the Fourier transform of 
iS(r^) is well valid while in high energy regime i.e. at 
small-a; and all the evolution characteristics for P{k±) 
exclusively contained in iS(rj_), so that we may write. 


and, 


P{k±,Y) = J d\x 5(rx,r) , 


dP{kx,Y) _ dS{k^,Y) 


dY 


= / dV_L e“ 


dY 
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( 8 ) 


and therefore from Eq. 
dq{Y) 1 


dY 


AttL 


dk\ 


,2^ '9‘5(r_L,U) 




dY 


= P\ 


[k_L,ri-] 


dS{rx.Y) 


dY 


( 9 ) 


Where for brevity in notations we have introduced P, 

= ^j dklkl J d?rx O(.10) 


It can also be shown that 

^[k±,r±] 

^[k±,r±] 


r]_ exp ( —^Lr\ 
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r_Lexp ( -^Lr± 


= 0 . 


( 11 ) 

( 12 ) 


Now, the non-linear evolution of the S'-matrix, iS(r_L = 
yl_—y'\)^ is governed by the Balitsky-Kovchegov equation 

w 


(BK) [2, [THSI, in large-A’c limit as, 

dS[yx,y'^-,Y) f ,2 (j/_L - 2 / 1 )^ 


dY 


d z±- 


27r2 J {yj_ - z±)'^(zj_ - y'^)^ 
[S{y±,y'±',Y) - S{y±, z±;Y)S{z±,y'j_,Y)] . 

(13) 























3 


Using Eq. @ and Eq. 0 with rj_ = ?/_L — 2/^, together 
with fact that the medium is transnationally invariant, 
evolution equation for q{Y), can now be written as, 


dq{Y) 

dY 


asNc 

2tt'^ 




(14) 


where A4{r±) is an integral over the daughter dipoles’ 
transverse coordinates. 


M(r^) = j d 


B^- 


(rx - B^YBl 
[S{r^,Y)-S{r^ - B^,Y)S{B^,Y)\ (15) 


The Balitsky-Kovchegov (BK) equation results from 
summing up long-lived soft gluon emissions, off almost 
onshell hard quark, where the gluon lifetimes usually 
much longer than the size of the target, which is taken 
to be a brick of size L here. The integral in Eq. ([6]) 
should therefore run from — oo to -foo, which means that 
the quark, scattering on the brick, is created at time, 
Ti = —OO, flies through semi-infinite empty space almost 
onshell while developing the gluon cascade giving the evo¬ 
lution of Eq. (USD, scatters on the brick, and flies off un¬ 
til Tf = -foo again developing the cascade m- Such a 
physical picture certainly does not apply in the context 
of jet quenching in heavy ion collisions where quark jets, 
are produced in the actual collision (at t = 0), with a 
high virtuality Q and quickly losses the virtuality by or¬ 
dered emission through (medium modified) DGLAP [H- 
evolution and have no time to develop a gluon cas¬ 
cade accordance with Eq. m before interacting with 
the plasma. In the present discussion we therefore pre¬ 
sume q as momentum broadening long after the quark 
leaves the brick when it again becomes almost onshell. 

In the event when U ~ 0 the jet quenching parameter 
qo and the saturation momentum Qso are approximately 
related as qo ^ AQ^^/L. In the rest of the article we 
will evaluate Eq. (HI inside or around saturation region 
which is the region of interest relevant in the context of 
transverse momentum broadening and energy loss phe¬ 
nomenology in relativistic heavy-ion collision studies. 


III. EVOLUTION OF q IN STOCHASTIC 
MULTIPLE SCATTERINGS REGIME 


Kinematic domains where in-medium interactions can 
be approximated by stochastic multiple scatterings, 
iS(r_L) appears to be a Gaussian in r± with variance 2/qL 

0. 


5[r_L] = exp 



(16) 


Above form of S{r±) in terms of qL and is limited, if 
not unique, in the sense that its Fourier transform leaves 
again a Gaussian form of P(fc_L) as given in Eq. JHl) and 
also return q once T acts over it, 





= q ■ 


(17) 


The imaginary part of the forward scattering ampli¬ 
tude A/'(rj_) in Glauber-Gribov-Mueller (GGM) multi¬ 
rescattering model also have a similar form, 

A/'(r_L) = l-exp|-^^T(6i)a;GAr ^a:, r^j , 

(18) 


where imaginary part of the forward scattering amplitude 
J\f is related with S'-matrix element, S{r±) as, 

Af{r±) = 1 - S{r±) . (19) 


As long as one is not deep inside the saturation region 
Eq. (fT8)) provides the correct form of saturation scale 
with the identification 

Qlib^) = ^Tib^)xGN ( 20 ) 

In order to calculate q at strong coupling via AdS/CFT 
correspondence, form of S as given in Eq. (HU) was even 
taken to define nonperturbative, quantum field theoretic 
definition of jet quenching parameter: q is the coefficient 
of Lr ^/4 in log S for small r± [l^ . In this paper we have 
assumed form of all S in Eq. (fT5|) as Gaussian in dipoles’ 
transverse size, with a variance of 2/qL, retains while it 
goes through high-energy evolution, with the caveat that 
deep inside the saturation region this may not true. We 
therefore take S'-matrices in Eq. (Ha as, 


5(rx,r) = exp 
5(r_L - B_l,F) = exp 
S{Bo_,Y) = exp 


m 

4 

m 

4 

m 

4 



U(r± - Bo_f 
LBl 


( 21 ) 

( 22 ) 

(23) 


Its important to recall that in this region of interest 
just inside the saturation region (and multiple stochas¬ 
tic scatterings works as well) S'-matrices are small, S <C 
1 and transverse width of daughter dipoles are large 
B_i_,Bj_ — r_L > 1/Qs- Once we replace Eq. (l^ - (l2^ 
in Eq. ([15]), M{rj_) becomes. 


A4(rj_) = exp 



d^Bi 


ir± - BA_fBl 


1 _ 


(24) 


However, the Gaussian assumptions of Eq. ([2T|) - Eq. 
(051) may not correct at low-a: both at small and large rj_. 
Small-a; evolution leads to anomalous dimension modify¬ 
ing the power of r_L from This means that q in Eq. 
(EH - Eq. (1231) is not a function of Y, but more appropri¬ 
ately is the initial condition of q{Y) at F = 0. Therefore, 
Eq. (l24l) should be understood as one step of BK evo¬ 
lution, where daughter dipoles interact through Glauber 
gluons are approximated by Gaussians in Eq. (1211) - Eq. 
(1551) . The subsequent evolution equations should then be 
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understood as a fairly crude approximation of the actual 
evolution. 

Now, around the saturation line, the dipole size rj_ ~ 
1/QsO, and B_\_, {B±_ — rj_) > \/Qs{Y) (however they are 
still well below \/Kqcd)- We will further assume that 
transverse size of the daughter dipole (~ inverse of the 
gluon’s transverse momentum) is much larger than the 
transverse size of the parent dipole i.e. B±, {B± — r±) ^ 
r±. With appropriate limits in the integral, Eq. ([M|) can 
then be written as, 


M{r±) 


exp 


m 


Lri 


r^/^QCD 


B/QI(y) 


Bi 


1 — e 


-^q(Y)LBl 


(25) 


Integral in Eq. (1251) is so converging that the upper limit 
may be taken at infinity instead of 1 /leaving the 
integral Idip as, 


Idip — 


dBl 

1 _ 

h/Ql{Y) Bl 



= Ql(X) p - exp 






_lq^\ 

2Ql{Y)) 


qL 




(26) 


where ‘Er’ is the Exponential Integral function. Using 
Eq. (l26l) one can now derive the evolution equation of 
the jet quenching parameter in stochastic multiple scat¬ 
terings regime just around saturation line as, 

dq{Y) asN, 2QUY) f q{Y)L \ 

dY TT U [ 2g2(y)J 

mL f q{Y)L V 
2QUY) A 2Q^s{Y))_' 

(27) 


Eq. (1271) can be written in terms of incomplete gamma 
fuction r(n, a), 


d In q(Y) 


= as 


ln2 + jE- 


^(l_e-2£)+r(0,2£:) 

oo 


(-£)* ' 

k{kl) 

(28) 


where. 


. 1 

4QUY) 


(29) 


with as = asNc/ir. Around f ^ 1, solution of Eq. (1^ 
can be approximately estimated as. 


/ , \ 0.9as 

q{x) oc f - j . (30) 


Jet quenching parameter g(x) therefore evolves with x, 
with an exponent 0.9 dg, which is just slightly less than 
that of xQ(x) where G(x) is the gluon distribution func¬ 
tion. 


A. Special Case: Double log enhancement 


For smaller dipole with r± ^ 1/Qs the imaginary part 
of forward scattering amplitude Af is small, and satura¬ 
tion and unitarity bound effects are not that very im¬ 
portant. Principle of color transparency ensures that the 
forward scattering amplitude JY approaches to zero (S- 
matrix, S, goes to one). This allows us to take that 
Af -C I for small dipole size even when small-a; evolution 
is included. For JY 1 we can linearise Eq. (HI as, 

In the limit of r± <C B±, Eq. (1511) can be further simpli¬ 
fied. 


M{r±) 



Lr\ exp 




B 


2 

±' 


(32) 


In a dipole system when a single real gluon emits, the 
gluon can be characterized by its energy, w, and trans¬ 
verse size, B±, that it makes with the (anti)-quark legs 
of the parent dipole having transverse size r±. When 
B± ^ r± the probability of emitting a gluon with given 
energy oj goes as r']_dB^/Bj_. Despite being non log¬ 
arithmic and converging, the integral become logarith¬ 
mically diverging once the single scattering in included 
which gives an additional B^ in the numerator This 
is evident from Eq. (|32D which together with the loga¬ 
rithm coming from energy w integral gives double log¬ 
arithmic contribution to the modification in transverse 
momentum broadening and so as for the jet quenching 
parameter. Substituting Eq. (HU Q-iid. Eq. (j32p in Eq. 
(fH|) gives the Y evolution of jet Quenching parameter, 

dq{Y) asN, fdBl 

Integral Eq. (1551) can be compare with Eq. (26) of Q 
and Eq. (4.35) of This double log enhancement, first 
derived in Q and subsequent later studies @-0, need to 
be supplemented with appropriate kinematic limits of the 
integral both for Bj_ and uj. 

In order to fix the kinematic limits for the ener^ w 
and B± here we follow the arguments made in to 
isolate the region of double log enhancement. Two con¬ 
dition have to be satisfied in order to achieve the double 
log enhancement. Firstly, the inverse transverse size of 
the daughter dipole 1/7?^ ("^ transverse momentum of 
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the emitted gluon) need to be just below saturation mo¬ 
mentum Qs but sufficient large that the multiple scat¬ 
tering is not important and one can linearise Eq. (USD 
with confidence. Secondly lifetime of the fluctuations (~ 
formation time of the emitted gluon) r ~ ujB\ have to 
be greater than nucleon size lo, and less than the length 
of the medium L. Inverse transverse size of the daugh¬ 
ter dipole is however well above the saturation scale Qs 
which ensure <C B\. 




1 1 
qlo ^ 


QCD 


Iq ^ ujB‘^ ^ L 


(34) 

(35) 


Eq. (ISSl) . together with the condition that Bj_ have to 
be sufficiently small so that one can linearise Eq. (USD, 
qLB\ < 1, implies, 


Bl 


< w < 



(36) 


As djdY = Lodjduj double logarithmic enhancement of 
jet quenching parameter in high energy would be, 

^9 = - / — / ^ 2 “ 

= f9(0)log^g, (37) 

where w integration have been performed before the B± 
integration. For an almost constant q{uj) we recovered in 
Eq. (IS71) the double log result (in the limit Qs —>• q^ 
first derived in Q and subsequent other studies [J@. 
This double log enhancement however diluted by multiple 
scattering effects 0 as evident from Eq. (E51) . 


introduced the concept of transverse deflection probabil¬ 
ity of a parton, that travels through strongly interacting 
medium, and derived high energy evolution equation for 
the jet quenching parameter in stochastic multiple scat¬ 
terings regime, which is the region of interest in the con¬ 
text of jet quenching phenomenology of the heavy-ion col¬ 
lider experiments. Balitsky-Kovchegov (BK) equation, as 
the evolution equation of the S-matrix, is used to derive 
high energy evolution equation for the jet quenching pa¬ 
rameter. We have shown that q{x) evolves with small 
X, with an exponent ^ 0.9 dg, which is just slightly less 
than that of xG{x) where G{x) is the gluon distribution 
function. The known result of double log enhancement 
emerges as a special case in the limit when the single 
scattering is only contributing. 

In this article Eq. (0 - Eq. (ITUl) constitute a complete 
set of equations needed to determine the energy depen¬ 
dence of q. To date there is no complete exact analytical 
solution of BK equation. Although there are several ap¬ 
proximate analytical solutions, along with some numer¬ 
ical solutions. One can therefore use the approximate 
analytical solutions, e.g, solution outside the saturation 
region which exhibits extended geometric scaling or solu¬ 
tion inside the saturation region, the Levin-Tuchin solu¬ 
tion [U , that exhibits complete geometric scaling [2^ to 
calculate the energy evolution of jet quenching parame¬ 
ter. As an obvious consequence of the scaling, one could 
end up finding energy-dependence is proportional to sat¬ 
uration momentum, Qs{Y), i.e., q(Y) ~ Q'^(Y). The re¬ 
sult could therefore modify at deep inside the saturation 
regions and also in the event when the LPM effect mod¬ 
ifies the spectrum strongly. Both the issues which may 
farther squeeze the evolution of jet quenching parameter, 
will be explore in a future effort. 


IV. SUMMARY AND OUTLOOK 

The are ongoing phenomenological efforts to extract 
values for the jet transport parameter q at various central 
heavy-ion collisions done at various energies for prevail¬ 
ing energy loss models. Here parameters for the medium 
properties are constrained by experimental data on the 
nuclear modification factor Raa [13 ■ Following a re¬ 
cent work by D’Eramo, Liu and Rajagopal [l| we have 
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